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SECTION  1.  INTRODUCTION 


The  Liquid  Wator  Content  Analysis  System  consists  of  three  major 
hardware  elements  and  a computer  program.  The  three  hardware  elements 
are  tlie  Scan  Converter,  the  Digital  Integrator  and  the  Inlerdata  7/32 
computer.  The  Scan  Converter  and  Digital  Integrator  were  developed  and 
delivered  on  previous  contracts.  Operation  of  the  Digital  Integrator  is 
explained  in  the  Equipment  Information  Report  for  the  Precision  Digital 
Video  Integrator,  Model  D-All-RS  dated  July  1974,  while  the  Scan  Converter 
IS  described  and  explained  in  the  Scan  Converter  for  Liquid  Water  Content 
Analyzer  Final  Report  dated  July  1975.  The  purpose  of  this  manual  is  to 
provide  information  essential  to  the  operation  of  the  system  through 
interaction  with  the  computer  and  the  computer  program. 


SECTION  2.  PKOCESSING  MODE 


As  shown  in  Figure  2-1,  the  computer  program  has  three  basic 
operating  modes:  lCECO.ttJ.3,  PEAYBACK  and  EIQUID  WAlEit.  CONI  EN  1 

ANALYSIS.  The  record  and  playback  modes  referred  to  in  this  respect 
are  record  only  and  playback  only.  The  liquid  water  content  analysis 
mode  is  the  primary  operating  mode  and  as  such  will  be  considered  in 
some  detail. 

In  the  paragraphs  that  follow  references  will  be  made  to  the  various 
constants,  variables,  tables  and  flags  used  by  the  computer  program.  In 
some  cases  it  may  be  necessary  to  refer  to  the  tables  in  Section  3 to 
understand  these  references.  All  teletype  messages,  allowed  responses 
and  subsequent  action  taken  by  the  program  are  presented  in  Section  4. 

A brief  description  of  all  subroutines  is  provided  in  Section  S,  while  more 
detailed  descriptions  of  all  major  subroutines  are  included  in  the 
Appendices. 

2.  1 RECOKD  Mode 

The  record  only  mode  is  used  to  record  the  raw  integrator  digital 
output  together  with  the  radar  ancillary  data  and  time  from  the  scan  con- 
verter. Data  is  recorded  in  the  form  of  table  VBUI'.  The  radar  calibration 
constant  associated  with  the  data  to  be  recorded  is  requested  on  the  tele- 
type before  the  recording  is  started.  The  teletype  is  armed  during  the 
recording  process  in  order  that  the  recording  can  be  stopped  by  typing 
HALT.  A tape  recorder  malfunction  or  end-of-tape  will  also  terminate 
the  recording. 

2.  2 PLAYBACK  Mode 

The  playback  only  mode  is  used  to  playback  data  recorded  in  either 
the  record  only  mode  or  the  liquid  water  content  analysis  mode  for  display 
on  the  scan  converter.  The  time  of  the  recording  and  the  integration 
parameters  used  during  the  recording  are  output  to  the  teletype  from  data 


Figure  2-1.  System  Flow  Diagram 


in  the  first  recorded  block.  1 he  teletype  is  armed  during  the  playback 
process  in  order  to  permit  termination  of  the  process  by  typing  HALT. 

A tape  recorder  malfunction  or  end-of-file  will  also  terminate  the  playback 
process. 

2.  3 LIQUm  WATEK  CONTENT  ANALYSIS  Mode 

Liquid  water  content  analysis  can  be  performed  on  either  real-time 
data  or  previously  recorded  data.  If  real-time  data  is  to  be  used,  the 
operator  can  choose  to  record  the  data  also.  After  the  input  mode  is 
established,  the  program  goes  into  a holding  loop  to  wait  for  teletype  input 
direction.  While  waiting  data  is  recorded,  played-back  or  ignored  at  a 
speed  determined  by  the  integrator.  Each  of  the  permissible  teletype 
commands  and  the  action  that  follows  them  are  summarized  in  the  following 
paragraphs. 

2.3.1  EDIT 

Typing  EDIT  will  permit  the  operator  to  change  any  of  the  acquisition 
or  analysis  parameters.  These  parameters  arc  divided  into  two  groups; 
the  operator  must  select  which  parameter  set  he  wishes  to  edit;  ACQ  or 
ANA. 

2.  3.  1.  1 Acquisition 

If  the  operator  selects  ACQ,  the  acquisition  parameters  are  re- 
quested along  with  their  format  in  the  following  order:  DIS,  PHIT , HCM, 

HM,  HO,  HE.  From  these  values  and  the  integration  parameters  the 
following  values  are  computed:  CPHIT,  SPHIT,  COPHIT,  HMH,  DH, 

KCK. 

The  program  then  builds  a table,  RLUT,  which  is  unique  to  this 
geometry.  Addressed  by  elevation  angle  PHI,  this  table  includes  the  range 
cell  number  R(PHI)  closest  to  the  trajectory  for  every  possible  PHI, 

Since  we  have  12-bit  angle  coding  but  elevation  is  restricted  to  one  quadrant, 
the  table  must  have  102‘1  entries.  Only  a small  fraction  of  the  entries  in 
KLUT  will  actually  be  used,  but  we  need  them  all  because  the  actual  values 
of  PHI  at  each  integrator  output  (dump)  are  not  known  apriori  since  the 
antenna  scan  is  not  synchronous  with  the  integrator  timing. 


(1) 


The  calculations  to  fill  HLUT.  described  by 

R(PHI)^DIi>/(COS  (PHI)  - SIN  (PHI)  COT  (PHIT)  ) 
involve  use  of  the  circular  COKDIC  subroutine  described  in  the  Appendix 
for  determination  of  the  transcendental  functions.  The  COl  (PHIl)  term 
needs  to  be  calculated  only  once,  while  the  SlN(PHl)  and  COS  (PHI)  terms 
must  be  obtained  for  each  PHI  in  the  Table.  The  COKDLC  routine  outputs 
both  SIN(PHI)  and  COS(PHI)  simultaneously. 

The  values  of  ■il.(PHI)  in  ilpHUT  are  scaled  in  teims  of  a ranpe  cell 
number  between  DIS  and  lOZi,  so  that  10-bit  words  are  required.  If  memory 
usage  IS  a problem,  schemes  for  reducing  this  word  length  to  8-bits  can 
be  envisioned  since  H(PHI)  must  be  monotonic.  lor  the  present,  Ib-bit 
half  words  will  be  used.  The  extra  bits  could  be  used  as  flags  to  show  where 
the  range  goes  beyond  1024  or  within  the  clutter  range  -KC. 

During  this  initialization  phase,  the  altitude  for  each  PHI  could  have 
been  easily  calculated.  However;  the  resulting  increased  data  handling  re- 
quirements in  the  acquisition  phase  to  follow  would  have  lengthened  the 
critical  execution  time  of  that  phase. 

Program  control  returns  to  the  holding  loop  to  await  furtlier  teletype 
input  instructions. 

2.3.  1.2  ANAly  sis 

If  the  operator  selects  ANA,  the  analysis  parameters  are  requested 
along  with  their  format  in  the  following  order:  DH,  HIl,  H12,  Z21NP,  HI3, 

Z3INP,  Tl,  T2,  T3,  SI,  S2,  S3,  S4,  Ql,  U2,  U3,  04,  DISNUM.  Acquisition 
parameters  are  then  recalculated  and  control  returned  to  the  holding  loop. 

2.  3.  2 PARAmeter  Display 

If  the  operator  selects  PARA,  the  program  will  print  all  acquisition 
and  analysis  parameters  on  the  line  printer.  In  this  manner  the  operator 
can  review  the  parameters  and  determined  which,  if  any,  of  them  need  to 
be  changed  by  EDITing.  ■' 

After  the  parameters  have  been  printed,  program  control  is  returned 
to  the  holding  loop  to  await  further  teletype  input  instructions. 
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Z.  3.  3 ACQuisition  Phase 

li  the  operator  is  satisfied  witli  the  parameters  that  have  been 
established  in  the  initialization  phase,  he  instructs  the  processor  (through 
the  TTY)  to  enter  the  acquisition  phase  whereupon  his  interaction  with  the 
machine  ceases  until  this  phase  has  been  completed. 

As  shown  in  the  flow  diagram  of  Figure  2-Z,  the  acquisition  phase 
begins  when  the  elevation  scan  reaches  one  or  the  other  of  its  turnaround 
limits.  A table  called  VBUF  is  periodically  being  refreslied  with  new  video 
data  for  each  of  10Z4  range  cells,  and  certain  corresponding  ancillary  data 
of  which  only  PHI  is  of  concern  here.  VBUh'  is  filled  via  DMA  transfers  from 
either  the  integrator  (real-time)  or  mag  tape.  The  contents  of  VBUF  at  a 
given  time  corresponds  to  one  of  the  sectors  shown  in  Figure  2-3.  The 
algorithm  within  the  loop  uses  data  from  KLUT  to  define  variable  length 
block  transfers  from  VBUF  to  anotlier  table,  VTAB,  For  ranges  below  KCi-C 
V(KCK)  is  entered  into  VTAB  instead  of  V(K),  When  the  antenna  reaches  the 
second  turnaround,  the  operator  is  instructed  that  the  acquisition  phase  has 
been  completed,  VTAB  now  contains  data  corresponding  to  the  shaded  area 
of  Figure  2-2;  these  data  are  inputs  for  the  analysis  phase  to  follow. 

Because  the  integrator  could  be  set  up  for  as  few  as  6-4  sweeps 
integrated  at  a PKT  of  9b0  sec.  , the  time  available  for  one  pass  through 
the  loop  in  Figure  2-2  is  bh.bb  milliseconds.  For  this  reason,  every 
effort  has  been  made  to  minimize  processing  within  the  loop;  even  at  the 
expense  of  increased  complexity  and  time  in  the  two  other  phases. 

I he  radar  data  associated  with  the  maximum  elevation  angle  greater 
than  PHIMAX  reached  during  the  acquisition  phase  is  saved  in  a table 
(KBHM)  to  lie  used  as  data  for  those  points  along  the  trajectory  above  ITM 
or  beyond  the  maximum  radar  range. 

When  acquisition  has  been  completed,  control  is  either  returned  to 
the  holding  loop  to  await  teletype  input  instruction  or,  if  not  in  the  real- 
time recording  mode,  control  is  passed  directing  to  the  ANALYSIS  phase 
software. 
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iguru  2-2.  Data  Acquisition 


Line  of  constant  PHI  where 
all  1024  range  cells  having 
been  integrated  since  the 
last  line  are  available. 


Launch  Site 


h'igure  2-1.  LWCA  Data  Acnuisition  Geometry.  Integrated  Log 
Video  for  the  Shaded  Area  is  Lntered  Into  VTAB. 
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l.i.4  ANALvsis 


Z.  3.4.  1 C (.nil putat ion 

i ho  analysis  phasi.’  Si.)it\varo  oompulos  the  Ik^uicI  \,vator  oonlont  density 
as  a function  (.>1  altitude  b.ised  on  the  values  of  radar  video  stored  in  the 
tables  V' 1 A B <ind  KBllM  and  the  analysis  pii'ameters. 

The  program  first  computes  the  altitude  extremes  for  which  there 
are  data  and  then  avi>rages  the  video  over  the  altitude  resolution  steps  (1)11) 
between  these  extremes.  1 ransition  altitudes  are  then  computed  in  units  of 
altitude  steps  and,  for  112  and  111,  as  a function  of  Z2INF  and  Z5INP, 
respectively.  The  coefficients  (U)  and  exponents  (fa)  for  eai.h  altitude  region 
are  next  interpolated  across  each  of  the  transition  regions  1 1.  12,  and  I 1. 

Values  ol  the  logarithm  ol  the  liquid  water  content  densitv  and 
liquid  water  content  density  are  computed  in  the  subroutine  l.WC  tor  each 
altitude  step  and  saved  in  tables  LW  1 AB  and  WlAB,  re  spect  i viT\ . 

The  calculations  used  to  deterniine  the  entries  for  t.tbli-s  1,'iV  lAB 
and  Wl  AB  are  presented  in  the  appendices  titled  "AnaUtic  Determination 
of  Liquid  Water  Content".  "Range  of  Lujuid  Water  C \.ntent".  aiui  LUt  faubroutine 

2.  3.  4.  2 Out  put 

When  the  liquid  wiator  content  calculations  have  been  completed, 
control  IS  transferred  to  a series  of  program  routines  for  outputting  the 
analysis  data.  1 he  first  of  these  routines,  LWCPLT,  outputs  the  data 
from  LW  1 AB  to  the  scan  converter  display  designated  by  DIfaNL'M.  (A 
programming  description  for  display  interface  programs  is  included  in  the 
appendix  titled  "Display  Data  Port  Programming".)  After  the  plot  of 
liquid  water  content  versus  altitude  has  been  presented  on  the  scan  convertor 
display,  the  program  outputs  all  acquisition  and  analysis  parameters  to  the 
line  printer  and  then  prints  a table  of  averaged  radar  video  and  calculated 
liquid  water  content  density  (WlAB)  on  the  line  printer. 

Program  control  is  then  transferred  back  to  the  holding  loop  to 
await  further  teletype  input  instructions. 


SECTION  3.  CONSTANTS,  VAEJABLES,  FLAGS  AND  TABLES 


3.  1 Constants  and  Vafiables 

Table  3-1  lists  all  the  significant  program  constants  and  variables. 

For  each  nanied  constant/ variable  there  is  a brief  description,  a range  of 
values,  the  unit  of  the  number,  the  source  and  the  scaling  of  the  constant 
or  variable. 

3.  2 Flags 

lable  3-2  on  page  3-11  lists  the  significant  program  flags  together 
with  their  use,  size,  number  of  states  and  state  definitions. 

3.3  Tables 

Program  data  tables  are  described  in  Figures  3-1  through  3-8. 
Applicable  notes  are  included  with  each  figure. 
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Table  3-1.  Coiistanls  imd  Variables 
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fable  3-1.  (Continuecli 


Multiplier 

Name  Description  R-m^e  W ithuut  Scaline  Units  Source:  Format  Belore  Storage  Value  of  LSB 


Use  Size  No,  of  Stales  State  Definition 


states;  there  a 
several  status 
cond  itions) 


Address  VBUF 

bi  f I 


bit 

31^ 


Address  - 
VBUF  ^-32 


VBUK  is  the  table  which  contains  the  integrated  video  data  and  its 
associated  radar  ancillary  data.  If  real-time  data  is  being  used,  the  video 
and  integrator  parameters  come  from  the  integrator  and  the  date,  time  and 
antenna  position  come  from  the  ancillary  data  port.  If  tape-recorded  data 
is  being  used,  all  data  comes  from  magnetic  tape.  The  tape  record  is 
formatted  just  as  VnUK.  Unit  s - seconds,  a four-bit  BCD  value  is  stored  as 
the  least  significant  four  bits  of  the  second  byte.  This  permits  determination 
of  time  to  the  nearest  ten  seconds  by  examination  of  only  the  second  half 
word  of  the  table. 


Address  - RLUT 


bi  t \ 


^ Address  - RLUT  + 204b 

Figure  3-Z,  R^inge  Look-Up  Table  (RLUT) 

RLUT  is  a look-up  table  which  identifies  the  range-cell  number 
which  lies  on  the  rocket  trajectory,  as  a function  of  radar  antenna 
elevation  angle.  There  are  1024  possible  antenna  elevation  angles,  PHI, 
and  a maximum  of  1024  radar  range  cells.  Each  R(PHI)  is  therefore 
a ten-bit  value. 

R(PHI)  = DR/(cos(PHl)  - sin(PHl)  . cot(PHIT)  ) 


Figure  3-3.  iiange  Cells  Beyond  HM 


RBHM  is  a table  containing  the 
for  all  range  cells  associated  with  the 
maximum  antenna  elevation  angle. 


integrator  (ta[ 
integration  pel 
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bit  ^ Address  = ADTABI 


bit 

31 


- J 

SEC 

TIME  1 

1 THETA 

PHI  ! 

Address  = ADTABI  6 


bit  Address  = AD  TAB2 


bit 

31 


L 


t 


DAY 

SEC 

TIME 

theta 

, PHI 

Address  = ADTAB2  t-  6 


figure  3-4.  Ancillary  Data  Tables 


The  ancillary  data  tables  contain  the  radar  ancillary  data 
associated  with  the  integration  periods  occurring  at  the  maximum 
and  minimum  elevation  angles.  Data  is  transferred  to  ADTABI  at 
the  first  antenna  turn-around  point  and  to  ADT  AB2  the  second  time 
the  antenna  scan  direction  changes.  Units- seconds,  a four-bit  BCD 
value,  IS  stored  as  the  least  significant  four  bits  of  the  second  byte 
in  each  case. 
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bit  Address  - VTAB  bit 


Figure  3-5.  Video  Table  (VTAB) 


VTAB  is  a table  of  selected  values  of  integrated  radar  video.  The 
selected  values  are  stored  in  VTAB  in  an  altitude  ordered  and  indexed 
manner.  The  index  for  storage  into  VTAB  is  computed  as  the  range-cell 
number  times  the  sine  of  the  elevation  angle  times  eight.  (This  index  is 
related  to  altitude  as  follows:  Altitude  - Index  * Range  cell  size/ 8).  The 

range  cell  values  chosen  for  placement  in  V'l'AB  are  chosen  on  the  basis 
of  range-cell  intersection  with  or  proximity  to  the  rocket  trajectory.  The 
determination  of  the  appropriately  ranged  video  was  made  by  reference  to 
KLUT  after  each  integration  period. 

The  range  values  chosen  for  the  beginning  table  entries  are  the  video 
values  at  the  ground  clutter  range  (KCK)  if  the  radar  range  to  the  trajectory 
intersection  is  less  than  HCK. 

The  maximum  index  into  V’J  AB  is  limited  to  159V.  This  limit 
corresponds  to  a minimum  altitude  of  15  kilometers. 


Figure  3-6.  Averaged  Video  Table  (VIDAVG) 


VLDAVG  contains  averaged  values  of  radar  video  for  each  of  the 
altitude  steps.  The  value  of  each  table  entry  is  between  0 and  511 
(9-bits)  and  therefore  each  table  entry  requires  Z bytes. 

The  number  of  table  entries  is  a function  of  the  altitude  quanti- 
zation, DH,  and  the  altitude  extremes  for  which  LWC  will  be  calculated, 
HO  and  HF.  The  first  entry  would  be  the  value  for  the  first  altitude 
quantization  step  above  HO,  the  second  entry  will  be  for  the  second 
step,  etc.  The  maximum  number  of  entries  will  be  Z56. 
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Figure  3-7.  Logaritlim  of  Liquid  Water  Content  Density  (LWTAB) 

LWTAB  contains  scaled  values  of  log-base-10  of  the  calculated 
liquid  water  content  density  (W).  Value  of  LW(H)  is  between  0 and  255 
to  permit  convenient  output  to  the  display. 

The  number  of  table  entries  is  a function  of  the  altitude  quanti- 
zation, DH,  and  the  altitude  extremes  for  wliich  W was  calculated,  HO 
and  HE.  The  first  entry  would  be  the  value  for  the  first  altitude  quan- 
tization step  above  HO,  the  second  entry  will  be  for  the  second  step,  etc. 
The  maximum  number  of  entries  will  be  256. 
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SECTION  4.  TEEETYPE  MESSAGES 


All  of  the  toletypo  messages  and  their  purposes  and  appropriate 
operator  responses  are  listed  in  Table  4-1.  'I  he  messages  are  ordered 
and  grouped  as  they  might  appear  on  the  teletype  machine.  Some  message 
can  appear  only  with  others  and  these  messages  are  grouped  between 
double  horizontal  lines  in  the  table,  lor  example,  all  of  the  acquisition 
parameters  become  available  for  editing  together  and  these  messages  are 
therefore  grouped  between  double  lines. 


NNNNN  METER  Permits  entry  of  nominal  third  traniition  Enter  transition  altitude  in  meters 

altitude  CR  will  force  use  of  previous  value 


SECTION  5.  PKOGiOVM  MODULES 


'I  he  complete  computer  progriim  consists  of  eight  (b)  program 
modules.  Each  of  these  modules  is  assembled  as  a stand-alone  program 
witli  references  to  common  constants,  variables,  subroutines  and  entry- 
points  in  other  program  modules  established  by  a series  of  equate  state- 
ments at  the  beginning  of  the  program  module.  '1  he  functions  performed 
by  each  of  the  individual  modules  are  described  in  the  following  paragraphs. 

Table  5-1  lists  all  important  subroutines  and  subprograms  together  with 

a brief  description  and  access  information.  1 able  5-1  begins  on  page  5-5.  ' 

I 

5.  1 Tables  and  Constants  (TBL) 

The  first  program  module  establishes  the  memory  location  and,  j 

if  applicable,  the  initial  value  of  all  program  constants  and  variables  which 
are  shared  among  tlie  other  seven  program  modules.  In  addition,  two  very 
short  subroutines  are  included  in  this  program  modules:  TTYEiiKOrt  is  a 

general  purpose  routine  which  simply  prints  EKKOri  on  the  teletype  machine 
and  then  returns;  IIYKCC  is  called  to  input  the  radar  calibration  constant 
from  the  teletype,  convert  the  input  to  binary  and  store  the  converted  and 
scaled  constant  for  future  use. 

5.  Z Executive  (EXC) 

The  executive  program  module  contains  the  routines  for:  1)  initial- 

ization of  program  flags  and  machine  interrupts;  Z)  determination  of 
processing  mode;  3)  the  KECOrUD  only  mode;  4)  the  PLAYBACK  only 
mode;  5)  console  display  of  the  radar  antenna  azimuth  and  elevation;  b) 
the  licjuid  water  content  holding  position  while  waiting  for  specific  analysis 
instructions;  7)  the  liquid  water  content  data  acquisition  phase;  8)  magnetic 
tape  command  generation  and  tape  recorder  status  processing;  9)  integration 
parameter  acquisition,  de-coding,  storage  and  output  to  the  teletype;  10) 
conversion  of  coded  time-of-day  to  ASCII  and  its  output  to  the  teletype; 
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li)  setting-up  the  scan  direction  iiags  for  data  acquisition;  12)  arming 
the  teletype  lc>r  input  command  while  performing  other  processing  tasks; 
and  13)  processing  radar  input/ output  interface  errors. 

The  entry  points  from  other  program  modules  are  INTl'l  for  re- 
initialisation; ACQCHtCK,  the  LWC  analysis  mode  holding  looping;  and 
lIMlfTTY  for  converting  the  coded  time-of-day  to  ASCII  and  outputting  it 
to  the  teletype. 

5.  3 £.dit  (EDT) 

The  edit  program  module  contains  all  the  processing  associated  with 
the  input  of  the  acquisition  and  analysis  parameters  used  tor  liquid  water 
content  analysis.  There  are  three  program  entry  points;  ANAEDIT  is  the 
entry  point  for  editing  the  analysis  parameters,  EW2  is  the  entry  point  for 
editing  the  acquisition  parameters  and  LW2D2  is  the  entry  point  for  rebuilding 
the  tables  and  re-computing  variables  that  depend  on  launch  trajectory 
parameter  s . 

5.  4 Analysis  (ANL) 

The  analysis  program  module  calculates  the  liquid  water  content 
density  along  the  trajectory  on  the  basis  of  the  radar  data  saved  in  VTAB 
during  data  acquisition  and  the  analysis  parameters.  There  is  a single- 
entry point  (ANALYSIS)  and  the  exit  is  to  the  output  program  module 
(ANAPLT)  whicli  outputs  the  calculated  analysis  values. 

5.  5 Subroutines  (SUB) 

Seven  major  program  subroutines  which  are  called  from  the  other 
program  modules  are  contained  in  the  subroutine  program  module. 

BCDASC  is  a subroutine  used  to  convert  binary-coded  decimal  (BCD) 
numbers  to  their  equivalent  ASCII  code  so  they  can  be  output  to  the  teletype 
or  line  printer. 

The  subroutine  ASCBIN  converts  a string  of  ASCII  numbers  to  a 
single  binary  (full  word)  equivalent.  The  binary  equivalent  is  scaled  up 
(multiplied)  by  a factor  equal  to  ten  to  the  number  of  digits  to  the  right  of 
a decimal  point. 
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Subroutine  BINASC  converts  a binary  lull  word  input  to  a string  oi 
ASCII  numbers.  If  the  input  is  negative,  the  most  significant  ASCII  character 

IS  a minus  sign;  if  the  input  is  positive,  the  most  significant  ASCU  character 
is  the  most  significant  digit. 

The  subroutines  CKOT  and  CVEC  are  circular  COKDIC  subroutines 
used  to  compute  trigonometric  functions  and  to  do  vectoring.  See  Appendix  A. 

The  teletype  driver,  TTYIO,  is  also  included  in  the  subroutine 
package.  This  subroutine  handles  all  teletype  input  and  output  on  an 
interrupt  serviced  basis. 

The  subroutine  FINISH  is  used  to  end  the  recording  process,  and 
write  an  end-of-file  (EOF)  on  the  magnetic  tape.  Control  is  return  to  the 
initiali^-ation  routine,  INITl . m the  executive  module. 

5 • b Magnetic  T ape  Driver  (T A P) 

The  magnetic  tape  driver  program  module  is  devoted  exclusively  to 
handling  input  and  output  for  the  magnetic  tape  recorder.  There  are  unique 
driver  comtnands  for  output,  input,  write  EOF,  rewind,  search  for  EOF, 
and  backspace.  Since  the  magnetic  tape  recorder  is  a 7-track  machine 
and  the  computer  memory  is  8-bit  oriented,  the  data  to  and  from  the  tape 
recorder  must  be  unpacked  and  packed  to  achieve  maximum  tape  utilization 
and  to  insure  no  loss  of  data;  the  packing  and  unpacking  routines  are 
included  as  part  of  the  magnetic  tape  driver, 

5 . 7 Kida r I/mj rivers  (HD K ) 

The  radar  I/O  driver  program  module  contains  the  subroutines  for 
the  input  and  output  of  radar  video  and  ancillary  data.  Video  data  are 
passed  to  and  from  the  computer  across  an  interface  with  the  integrator. 
Ancillary  data  (date,  time,  antenna  azimuth  and  elevation  angles)  are 
exchanged  through  an  interface  directly  with  the  scan  c onverter,.  refresh 
memory  system.  Video  data  input  and  output  are  handled  by  subroutines 
VIDINP  and  VIDOUT.  respectively;  ancillary  data  input  and  output  are 

handled  by  ANCINP  and  ANCOUT,  respectively,  bee  Appendix  G for  more 
details. 
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5.  8 Output  (OUT) 


Special  output  routines  for  liquid  water  content  data  are  included  as 
part  of  the  output  program  module.  The  primary  entry  point  is  to  ANAPLT 
This  entry  is  used  when  the  results  of  liquid  water  content  analysis  are  to 
be  plotted  and  printed  in  tabular  form.  The  routine  for  printing  all  of  the 
system  acquisition  and  analysis  parameters  is  also  part  of  this  program 
module.  Tliis  routine  (DIPAKM)  is  called  by  entry  to  .ANAPLT  but  niay  be 
called  directly  so  that  the  parameters  can  be  inspected  before  acquisition, 
analysis  or  editing. 
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rile  CORDIC  technique  ‘ ‘ for  computing  transcendental  functions 

has  been  chosen  for  implementation  of  certain  subroutines  required  in  the 
1-^^  C initialization  and  analysis  phase  software.  The  method  is  appealing 
because,  by  using  simple  add,  subtract,  shift,  and  com])are  operations 
which  are  available  as  computer  instructions,  the  look-up  table  require- 
ments are  reduced  dramatically  over  a pure  look-up  table  approach. 
Compared  to  a series  e.xpansion  method,  CORDIC  has  a speed  advantage^' 
and  it  calculates  two  functions  simultaneously.  Details  of  the  algorithms 
are  well-covered  in  the  r eferences  ^ ‘ and  will  not  be  repeated  here. 

The  CORDIC  algorithms  involve  iterative  operations  among  three 
registers  X,  and  X in  a process  similar  to  that  employed  in  successive 
approximation  A/D  converters.  The  schematic  representations  of  figure  1 
show  the  functions  which  can  be  calculated  by  means  of  the  two  circular 
CORDIC  subroutines,  named  CROT  and  CVKC.  In  the  rotation  mode, 

(CROT),  operations  are  performed  such  that  the  /,  (angle)  register  is  driven 
toward  zero,  while  in  the  vectoring  mode  (CVEC),  the  Y register  is  driven 
toward  zero.  The  constant  Kj  («=  1.  f>5),  sometimes  called  the  stretch  factor. 
IS  a function  only  of  the  number  of  iterations.  The  first  two  examples  of 
applications  which  appear  in  I'igurc  1 involve  straightforward  evaluation  of 
trigonometric  functions.  Example  3 demonstrates  how  the  multiple 
input/output  feature  can  be  used  to  advantage  in  evaluating  a more  complex 
trigonometric  expression.  Note  that  after  step  b,  register  X contains  the 
proper  number  for  the  input  at  step  d.  Similarly,  after  step  d,  register 
X contains  the  required  input  for  step  f.  In  cases  such  as  these,  the  numbers 
can  simply  be  left  in  the  registers  even  while  other  operations,  such  as  the 
addition  in  step  c,  are  executed. 
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The  scaling  for  sine,  cosine  and  angle  quantities  in  the  LWC 
software  has  been  established  as  described  in  reference  5.  The  sine 
and  cosine  are  scaled  alike  and  contain  15  significant  bits  such  that  a 
real-world  value  of  unity  corresponds  to  3Z7  68  (y'8000')  in  the  Interdata 
7 / 3Z  computer.  The  angle  is  scaled  so  that  90  degrees  is  represented  as 
10Z4  (y'400').  The  consensus  of  references  J thru  4 is  that  for  15-bit 
accuracy,  19-bit  registers  and  16  iterations  are  required.  The  entire 
3Z  bits  of  the  general  purpose  registers  could  have  been  used,  but  in 
order  to  ensure  at  least  19  bits  for  accuracy  while  allowing  capability  for 
numbers  much  larger  than  those  scaled  like  sines  and  cosines,  the 
following  arrangement  has  been  adopted.  The  X and  Y registers  are 
shifted  left  5 places  and  the  Z is  shifted  10  places  before  the  algorithm  is 
executed.  After  completion,  the  regi.sters  are  shifted  back  a like  number 
of  places.  Considering  the  stretch  factor,  inputs  to  X and  Y should  be 
kept  to  numbers  within  + to  avoid  overflowing  the  registers.  Inputs 
to  X and  Y can  be  scaled  arbitrarily  (what  goes  in  comes  out  multiplied 
by  K|),  but  numbers  which  are  smaller  than  those  scaled  like  sines  and 
cosines  result  in  increased  algorithmic  errors. 

Inputs  to  the  Z register  must  be  scaled  pro[>erly  since  the  look-up) 
table  for  tan  ^Z  * used  in  CROT  and  CVEC  has  predetermined  scaling. 

The  algorithms  will  work  for  values  of  Z over  a range  of  about  + Y*470' 

{+  99.  8°).  Outside  of  this  range,  CROT  "limits"--it  outputs  the  numbers 
corresponding  to  the  nearest  extreme.  To  prevent  overflow  of  the  Z 
register,  keep  inputs  within  + K 

The  CAL  assembly  listing  for  CROT  and  CVEC  has  been  included 
as  Figure  Z.  It  would  have  been  possible  to  code  CROT  and  CVEC  as  one 
program  with  its  mode  controlled  by  a flag  since  they  differ  only  in  the 
test  on  Z or  Y,  but  to  conserve  execution  time,  it  was  decided  to  generate 
two  separate  subroutines.  They  do,  however,  share  the  arctangent  look-up 
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table  CALFAT.  Execution  time  has  been  estimated  at  588  microseconds 
for  either  routine.  Note  that  the  location  counter  save  register  to  be 
used  with  the  BAL  instruction  is  F,  while  the  I/O  registers  are  B,  C 
and  D.  In  addition,  registers  7 thru  A are  used  internally  and  should  be 
saved,  if  required,  before  GROT  or  CVEC  is  called. 

Both  subroutines  were  run  under  OS33/1dT  with  a simple  dri\er 
added  to  generate  1/Kj,  load  the  X,  Y and  Z registers,  and  execute  the 
BAL  instruction.  The  registers  were  e.xamined  through  the  hex  display 
console  after  each  of  a nvunber  of  tests.  GROT  and  CVTsC  performed  as 
e.xpected. 
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The  Liquid  Water  Content  Determination  System  being  developed 
under  contract  DNA  00  J -75 -C -0050  will  use  the  digitized  log-video  output 
from  the  Raytheon  Model  D-A11-R5  Precision  Digital  Video  Integrator 
(PDVI)  to  derive  liquid  water  content  densi'y  (VV)  as  a function  of  altitude. 
This  derivation  will  be  made  on  the  basis  of  the  following  relationship: 

Z = U . W ^ 

where 

Z = radar  reflectivity 

W = liquid  water  content  density 

Q = phase  related  constant 

S = phase  related  constant 

During  the  data  acquisition  phase  of  the  liquid  water  content 
determination,  a table,  VTAB,  will  have  been  established  and  will  contain 
values  of  the  log -video  output  for  points  along  the  rocket  trajectory  of 
interest.  Through  a process  of  coordinate  transformation, the  altitude, 
including  a correction  for  earth  curvature,  will  be  determined  for  each 
of  the  table  entries.  The  total  altitude  extent  for  table  entries  will  be 
divided  into  altitude  groups  of  a size,  DH,  be  specified  by  the  operator 
and  then  all  values  of  radar  video  contained  in  a given  group  wTll  be 
averaged  to  yield  a single  value  of  radar  video  which  is  representative 
for  that  altitude  segment. 

The  digitized  log-video  output  of  the  integrator  is  related  to  Z 
in  the  following  manner: 

1 0 l.o  g Z = V - R C C 

where  V = log-video  and 

RCC=  the  radar  calibration 

The  relationship  to  be  used  for  the  determination  of  W is 
therefore, 

V - RCC  r log  Q + .S  log  W 
therefore,  ,, 

r-Y-^l<CC.  . Q-| 

, -1  10 

W = log — 

where  as  previously  stated  S and  U are  constants  related  to  the  particular 
phase  of  the  water. 
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A total  of  four  water  phases  will  be  considered:  rain,  large 

snow,  small  snow,  and  ice.  The  altitude  at  which  these  phase  changes 
take  place  ha  e been  designated  as  follows: 


HI 

- 

rain  to  large  snow- 

HZ 

- 

large  snow  to  small  snow 

H3 

- 

small  snow  to  ice  . 

The  thicknesses  of  these  phase  change  regions,  centered  on  HI,  HZ, 
and  H3,  have  been  designated  TI,  TZ,  and  T3  respectively.  The  values 
of  Q and  S to  be  used  in  the  determination  of  W have  been  designated  as 
follows : 

Ql,  SI  for  rain 

QZ,  SZ  for  large  snow 

Q3,  S3  for  snvall  snow  and 

Q4,  S4  for  ice. 

The  values  for  Q and  S will  be  linearly  interpolated  across  each  of  the 
transition  regions,  Tl,  TZ  and  T3.  The  interpolated  values  of  U and  S 
will  thevv  be  used  in  the  determination  of  W. 

The  phase  transition  altitude,  HI,  will  be  determined  by 
observation  of  the  "bright-band"  on  the  radar  display  and  will  therefore 
be  input  by  the  operator. 

The  phase  transition  altitudes,  HZ  and  H3,  will  be  derived  as 
follows : 


1.  A minimum  altitude  at  which  the  phase  change  can  take  place 
will  be  determined  from  radiosonde  data  and  input  by  the  operator, 
HIZ,  HI3. 

Z.  A mininumi  value  of  radar  reflectivity  will  also  be 
determined  and  specified  by  the  operator  for  each  transition 
point,  ZZ,  Z3. 

3.  The  transition  altitude  HZ  shall  be  that  altitude  above 
HIZ  where  the  radar  reflectivity  is  less  than  ZZ. 

4.  The  transition  altitude  H3  shall  be  that  altitude  above  H3 
where  the  reflectivity  is  less  than  Z3. 

Table  1 lists  the  constants  required  for  the  determination  of 
liquid  water  content,  their  definition  and  their  input  source.  In  cases 
where  the  source  has  been  designated  as  both  TTY  and  Display,  the 
intent  is  to  indicate  either  of  these  devices  as  input  sources. 


IU4 


Unclassified 
JI1T;75:J9 
<iO  January  1975 
Fane  3 of  1 


Table  J 

LWC  Determination  Constants 


iiource 


Name 

Definition 

1 Derived 
1 

TTY 

Di  s play 

Integrator 

Z 

Radar  reflectiv-ity 

X 

W 

Liquid  Water  Content  Density 

X 

Ul,  QZ,  03,  04 

Fha.se  related  constant 

X 

SI.  S3,  S3,  S4 

Fhasc  related  constant 

X 

VTAB 

Range  ordered  table  of  radar 
video. 

X 

V 

Single  range  value  of  log -video 

X 

RCC 

Radar  calibration  constant 

1 

. j 

X 

HI 

Rain/ large  snow  transition 
altitude. 

1 

1 

X 

j 

X 

HZ 

Large  snow/small  snow 
transition  altitude. 

X 

1 

H3 

Small  snow/Ice  transition 
altitude. 

X 

HIZ,  HI3 

Minimiun  altitudes  for  HZ  | 

and  H3. 

1 1 

1 

1 

X ! 

X 

ZZ,  Z3 

Minimum  reflectivities  for 
HZ  and  H3. 

X 

1 i 

T1 

Thickness  of  the  rain/large 
snow  transition  region 
centered  on  H 1 . 

X 

1 

X 

1 

TZ 

Thickness  of  the  large 
snow/small  snow  transition 
region  centered  on  HZ. 

1 

1 ^ 

1 

i 

X 

i 

1 

T3 

Thickness  of  the  small 
snow/ice  transition  region 
centered  on  H3. 

1 

X i 

X 

1 

1 

DH 

Size  of  altitude  groups.  | 

1 

1 

HO 

Ixiwest  altituile  for  which  I 

I..WC  is  desired.  | 

! 

1 

i 

' 1 

1 

HE 

Hijihost  altitude  for  which  j 

LWC  id  desired.  | 

1 

i 

1 

1 

1 1 

j 

1 1 
! 
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Reference:  1.  JHT:75:19,  "Analytic  Determination  of  Liquid 

Water  Content",  dated  Re\-.  1--1/J6/75. 
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3.  J1IT:75;37,  "LWC  Cons  tants  / Variables/ Tab  less", 
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This  memorandum  establishes  the  range  of  the  variables  used  to 
determine  liquid  water  content  from  the  integrator  output  values  of  log- 
video.  From  these  values,  the  range  extent  of  W'  is  then  determined. 
RCC  shall  be  the  Radar  Calibration  Constant  such  that 

V - RCC  = 10  log  Z 

where  V is  the  eight-bit  integrator  output  (MSB  = 50  dB)  and  Z is  the 
radar  reflectivity. 

Example : 

68  dB  (1  0 1 0 1 1 10) 

48  dB 


V = 
RCC  = 


, -1  /68  - 48\ 


= 100 

From  JHT:75:19,  Revision  1, 
V - RCC 

W . - ' 


log 


10 


- log  Q 


where  W is  liquid  water  content  in  units  of  granis/m^. 

As  per  K.  M.  Glover  of  AFCRL  (LYW)  on  34  January  1075,,  the 
values  of  U and  S to  be  encountered  in  the  operational  envirotiment  are 
as  follows: 

100  ^ Q ^ 300,  000 

1.  63  i S s 54 

I’l'.i'u.t.Dr  T POE  3L.11K-:';CT  p-fL-l-lD 
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The  ritn.ue  of  V is 

0 ^ V < 100 
and  the  range  of  RCC  is 

0 S KCC  s 00 

The  analytical  range  of  W may  now  be  determined  as 
follow’ s : 

V - RCC 

.1  I 

VV  = log 
W is  maximum  when 
V - RCC 


10 


- log  Q 


S 


10 


log  O 


is  maximum. 


\V  = log 

max 


- 1 


'V  - RCC 
max  mm 

To 


loc  C 

min 


= log 


- 1 


'100-0  , 

. 1.  63 


m0‘.‘0  = 


W 

iTun 


log 


j - log  ZOO.  000 

1.  63 


4.7Z  X 10 


- 10 


The  range  of  VV  is  therefore, 


4.7Z  X lO"**'^  *•  VV  i MO'Hi'V 
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On  the  basis  of  practical  considerations,  this  range  can  be 
reduced  to 

10'^  s W ^ 100 

Such  a range,  5 orders  of  magnitude,  rules  out  the  possibility  of 
graphically  representing  W on  a linear  scale.  A logarithmic  scale  will 
be  used  for  plotting  the  \alues  of  W.  The  X-axis  of  the  display  will  be 
divided  into  five  decades  for  this  purpose  ranging  from  10'^  to  100. 

In  terms  of  software  this  will  mean  the  creation  of  a new  table, 
LVV'TAB,  which  will  contain  values  ol  log  VV  to  be  used  for  display  output. 

So  as  to  insure  that  all  arithmetic  manipulations  will  involve  only 
integers,  the  following  scaling  will  be  accomplished  in  the  course  of  the 
determination  of  log  W and  W: 

4.  TIve  quantity,  V RGG, — t-s  oonstrainod  to  be  g-feater  than 

tro 

Z.  The  quantity,  V - RCC,  will  be  multiplied  by  100  prior  to 
division  by  10;  i.  e.  , multiplied  by  10. 

3.  The  quantity  calculated  l)y  steps  1 and  i above  will  be 
increased  by  the  addition  of  three  (3)  times  the  stored  value 
of  S.  (The  stored  value  of  S is  100  times  the  real  world  value 
of  S.  See  reference  3.  ) 

4.  The  quantity,  log  Q,  will  always  be  positive  due  to  the 
permissible  range  of  G (400  - 200,  000)  and  will  be  multiplied 
by  100  prior  to  subtraction  from  the  quantity  established  in 
steps  1,  2 and  3. 

5.  The  value  thus  established  will  be  multiplied  by  the  quantity 
256/5  to  provide  for  display  scaling. 

6.  The  quantity  thus  established  the  steps  1 through  5 above 
will  be  divided  by  the  stored  value  of  S. 
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Tlicse  steps  will  yield 


100  • (V  - RCC)  ^ j 
10 


100 


10(1-  S 


Z56 

5 


(V  - RCC) 

To 


log  Q 


Z 56 
~5~ 


Had  step  five  been  omitted  (multiply  by  Z56/5),  the  result  v/ould  have 
been  the  addition  of  the  quantity  3 to  the  \ alue  of  log  W.  This  addition 
scales  the  range  of  log  W,  -3  to  tZ,  to  be  0 to  +5.  The  multiplier  of 
step  5 converts  the  log  W + 3 to  an  integer  between  0 and  Z56  in  antici- 
pation of  output  to  the  display. 

For  the  calculation  of  W from  log  W,  the  scaling  will  be  changed 
as  required  by  the  hyperbolic  cordic  function.  Special  scaling  of  the 
inputs  to  the  cordic  functions  used  to  calculate  logarithm  and  inverse- 
logarithm  will  be  required  as  described  in  A.  J.  .lagodnik's  mimiurandum, 
.AJJ-Z4. 
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Eiquid  Water  Content  Subroutine 
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LW  C Subroutine 


This  subroutine  solves  the  equation  for  liquid  water  content 


W 


-1  ^ 
log  — 


(1) 


where  V'  - RCC  (=  10  log  /.)  is  log  Video  corrected  for  the  radar  calibration 
constant  RCC,  Q is  the  phase- related  coefficient  and  S is  the  phase-related 
exponent  in 

Z = QW^  . (^) 

The  subroutine  makes  use  of  two  hyperbolic  CORDIC  algorithms  described 
in  AJJ-<d9  of  Appendix  id  HVEC  to  obtain  log  Q and  HROT  to  evaluate  the 
exponential  in  (1). 

The  inputs  are  left  in  the  following  registers  prior  to  calling  LWC: 

R3:  Average  Video  Z.  5bV  (From  LW  TAB),  0 < \’  < 99 

R4:  Exponent  100  S,  1.  63  < S < Z.  64 

R5;  Coefficient  Q,  100  <Q  < ZOOOOO 

R9:  Address  in  LWTAR 

The  subroutine  expects  to  find  the  radar  calibration  constant  at  \ RCC.  The 
outputs  remain  in  registers 

Rb;  (Log  W I Log  W scaled  for  LW'CPLT  (color  display) 

R7;  lOOOW’  ; W'  scaled  for  ANAPLT  (line  printer). 

When  LWC  execution  is  complete,  control  is  returned  to  the  address  in  R8. 
Although  other  registers  are  used,  their  contents  are  stored  and  replaced; 
only  the  contents  of  R6  and  R7  change  as  a result  of  calling  LWC.  In  addition 
to  the  register  outputs,  average  video  is  stored  in  the  V'IDAV'G  table,  indexed 
by  R9. 

The  input  O in  RS  is  tested  and  limited  to  the  range  mentioned  above, 

Z 3 Z4 

then  it  is  shifted  left  until  its  value  is  between  Z and  Z , with  the  number 

R 1 

of  shifts  required  (6  to  17)  being  contained  in  Rl.  The  number  in  R5,  Q Z 

becomes  the  variable  V in  example  4,  Figure  1 of  .AJ.T-Zd  in  Appendix  E 

Z4 

scaled  so  that  V = 1 corresponds  to  Z in  the  computer.  The  quantities  V 4 1 
and  V - 1 arc  loaded  into  the  X and  Y input  registers,  then  PVEC  is  executed. 
.Although  HVEC  and  HROT  were  written  as  subroutines,  they  arc  coded  in-line 
within  LWC  because  they  are  not  used  elsewhere.  When  HROT  has  ended, 

n-3 
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2 In  V remains  in  the  Z register  which  is  subsequently  multiplied  by 

200  log  e • 2““  so  that  the  n^ost  significant  )3art  of  the  result  ends  up  in  the 

Y register  as  100  log  e . In  V • 2^‘\  Next,  HI  is  similarly  multiplied  so 

that  R1  . 100  In  2 • log  e • resides  in  RO.  Finally,  RO  is  subtracted  from 

the  constant  Y'2D278D-1V  ( 1 ()0  log  e . In  2^'^.  is  added  t.;  Y so  that  its 

2 0 

\’alue  becomes  100  Log  Q • 2 

The  averaged  video  in  R3,  after  correction  for  the  similarly  scaled 

calibration  constant  RCC,  is  limited  to_^the  range  -77  < R3  < 286  and  multiplied 

by  a c.nistant  to  heccune  1 0(V  - RCC)  2^  . Next,  Y which  still  contains 

100  log  (J  • 2 is  subtracted  from  R3  and  the  result  is  divided  bv  1 00  • S in 

2 0 

R-4,  leaving  log  W • Z"  in  register  Z which  is  then  limited  to  the  range 


2^^<  Z<  2 


,20 


The  Z register  contents  are  then  duplicated  in  output 
register  Rn  and  a constant  added  to  develop  the  final  output  (log  VV  + 3)  2*^ 
to  be  used  in  LWCPLT. 

The  number  in  Z is  arithmetically  altered  by  constants  to  change  its 

base  and  rescale  it  so  that  (In  W - 6 In  2)  2^“^  ends  up  in  the  Z register.  The 
24 

constant  (In  2*2  ) is  added  to  Z and  the  number  of  additions  counted  in  R1 

until  Z is  greater  than  (1  - In  2)  Z^^.  Now  Z contains  (In  \V+  (R1  - f>)  In  2)  2^  \ 

\ is  zeroed  and  the  constant  2^"^/  (K  - 1)  is  loaded  into  Y so  that  the  exponential 
can  be  evaluated  using  HROT  as  in  example  3,  Figure  1 of  AJJ-29  in  .Yppendix  E. 
After  HROT.  the  addition  of  the  sinh  and  cosh  in  X,  and  a shift  left  of  six  bits, 


,R1 


,24 


that  register  contains  W . 2“'  . 2‘'"’.  Additional  shifting  by  R1  bits  in  the 

opposite  direction  yields  W • 2^  in  X which  is  rescaled  to  1000  • W and  left 
in  R7  as  the  output  for  ANAl^LT. 
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The  CORDIC  technique  ' for  computing  tranaendental  functions 
has  been  chosen  for  implementation  of  certain  subroutines  required  in  the 
LW C initialization  and  analysis  phase  softvvare.  Circular  CORDIC  sub- 
routines CROT  and  CV'EC,  us^d  in  trigonometric  operations  in  the  LW C 
softw'arei  have  been  described.  This  memo  considers  Hyperbolic 
CORDIC  subroutines  HR  OT  and  HVEC,  needed  in  the  LWC  analyzer  for 
evaluatioii  of  natural  logs  and  antilogs.  Hyperbolic  CORDIC  differs  from  the 
more-easily-visualized  circular  algorithms  in  that  the  vectors  involved  are 
constrained  by  a hyperbola  rather  than  a circle.  The  hyperbolic  CORDIC 
algorithms  share  the  same  advantages  over  alternative  techniques  as  the 
circular  algorithms:  1)  They  employ  add,  subtract,  shift  and  compare 
operations  available  as  rapidly-executable  computer  instructions,  2)  The 
amount  of  memory  required  for  stored  constants  is  minimal,  and  3)  Two 
functions  are  obtained  simultaneously. 

Authors  of  CORDIC  papers  routinely  point  out  how  the  circular 
algorithms  can  be  converted  to  hyperbolic  form  simply  by  changing  the 
stored  constants  and  decision  criteria.  In  actual  practice,  however, 
convergence  problems  require  special  attention  not  needed  in  the  case  of 
circular  algorithms.  The  only  treatment  of  these  problem.s  known  to 
this  author  is  to  be  found  in  Reference  4 which  also  suggests  a solution 
that  works.  TH  s comprehensive  paper-which  unifies  algorithms  elsewhere 
referred  to  as  circular,  linear,  and  hyperbolic  CORDIC-also  contains  val- 
uable information  as  to  convergence  domain  and  accuracy. 

The  CORDIC  algorithms  involve  iterative  operations  among  three 
reeistej-s  X.Y,  and  Z in  a process  similar  to  that  employed  in  successive 
approximation  A/D  converters.  Details  of  the  hyperbolic  algorithms  appear 
in  Refernce  4 and  will  not  be  repeated  here  except  to  note  that  repetition  of 
the  fourth  and  thirteenth  iterations  was  necessary  to  solve  the  previously 
mentioned  convergence  problem. 

The  schematic  representations  of  Figure  1 show  the  functions  which 
can  be  calculated  by  means  of  the  two  hyperbolic  CORDIC  subroutines, 
named  HRCTT  and  HVEC.  In  the  rotation  mode  (HROT),  operations  are 
performed  such  that  the  Z (hyperbolic  angle)  register  is  driven  toward 
zero,  while  in  the  vectoring  mode  (HVEC)  the  Y register  is  driven  toward 
zero.  The  constant  K|  is  a function  only  of  the  number  of  iterations  but 
must  take  into  account  the  two  which  are  repeated.  For  18  iterations,  K . has 
been  calculated  as  0,  828159361  on  the  9820A.  Note  that  this  number  is 
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about  half  of  the  constant  for  circular  CORDIC. 

The  first  t\\o  examples  which  appear  in  Figure  1 involve  straight  - 
forward  evaluation  of  hyperbolic  functions;  note  the  convergence  domain  lim- 
itations. It  is  interesting  that  X must  be  kept  non-negative  (observed  ex- 
erimentally).  If  the  number  in  the  X register  is  negative,  then  the  hyperbola 
shifts  to  the  other  half  plane  (see  Figure  1 of  Reference  4)  and  the  decision 
criteria  cause  rotations  in  exactly  the  wrong  directions,  so  that  the  Z register 
always  ends  up  at  one  or  the  other  of  its  saturation  limits.  Example  three 
illustrates  the  simplicity  of  exponential  evaluation  once  HROT  has  been  run. 
Finding  logarithms,  as  in  example  4,  is  a bit  more  tricky  and  makes  use 
of  an  identity  which  requires  addition  and  subtraction  before  HVEC  is  called. 
Note,  however,  that  only  a shift  by  one  bit  and  no  multiplication  to  correct  for 
K_|is  needed  since  it  cancels  out  in  register  Z.  A side  benefit  of  example 
4 is  that  it  also  obtains  after  a corrective  multiplication. 

HROT  and  HVEC  (Figure  2)  were  intended  to  have  the  same  15-bit 

precision  as  that  obtained  in  CROT  and  CVEC,  but  two  more  iterations  were 

added  (because  of  the  convergence  problem).  Thus  the  hyperbolic  routines 
employ  18  iterations,  and  at  least  19  -ignificant  bits  must  be  maintained  in 
the  three  I/O  registers.  Inputs  to  X and  Y can  be  scaled  arbitrarily  (what 
goes  in  comes  out  multiplied  by  K_j[),  but  numbers  with  fewer  than  19 
significant  bits  will  result  in  increased  algorithmic  errors. 

Inputs  to  the  Z register  must  be  scaled  properly  since  the  look-up 
table  HALFAT  has  prej-determined  scaling  where  a real-world  value  of 
unity  corresponds  to  2 or  Y'lOO  0000'  in  the  Ir.terdata  7/32  computer. 

HROT  and  HVEC  , icoded  separately  to  speed  execution  even  though 

they  differ  only  in  the  test  on  Z,  were  originally  written  as  separate  sub- 
routines sharing  the  table  HALFAT.  When  it  was  determined  that  each  would 
be  called  only  once  during  another  subroutine  (LWC),  they  were  both  in- 
cluded directly  as  part  of  that  subroutine.  Hence  the  missing  addresses  in 
Figure  2-the  CAL  assembly  listing  of  HVEC,  HROT,  and  HALFAT. 

The  fact  that  the  Fourth  and  thirteenth  iterations  need  to  be  re- 
peated meant  that  the  number  of  places  to  be  shifted  was  not  simply 
related  to  the  index  as  in  the  circular  algorithms.  Rather  than  slow 
execution  by  testing  I and  branching,  it  was  decided  to  include  a shift 
number  with  each  of  the  18  entries  in  hL^LFAT  which  had  sufficient  room 
available.  Thus,  each  full  word  in  HALFAT  (figure  2b)  contains  the  shift 
number  in  right,  most , byte,^and  the  most  significant  24  bits  of  the 
constant.  2^  tanh  ' (2  ^ ° ‘ in  the  remaining  bytes.  During  each 

iteration,  the  shift  number  is  put  into  SAR  by  means  of  a load  Byte  in- 
struction. After  the  shifts,  SAR  is  again  loaded  from  HALFAT  but  with 
the  full  wrord  this  time.  SAR  is  then  shifted  right  by  eight  bits  to  eliminate 
the  shift  humber  and  properly  scale  the  constant.  Even  though  the  first 
entry  in  HALFAT  is  a negative  number,  the  required  positive  quantity  in 
SAR  results  because  a logical  shift  right  instruction  was  used. 

Execution  time  for  either  hyperbolic  CORDIC  routine  has  been 
estimated  as  663  microseconds.  Test  calculations  of  sinh,  cosh,  and  Tanh 
at  various  points  scattered  over  the  convergence  domain  revealed  a peak 
error  corresponding  to  15-bit  accuracy  in  the  outputs. 
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GTO  AGAINC 
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Figure  2a.  H\’EC  Assembly  Listing  (Part  of  LWO 
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The  Display  Data  Interface  design  plan  contains  sections  entitled 
"Operation  of  the  LWCA  Control  Panel"  and  "I lardware/ Software  Interaction". 
The  purpose  of  this  memo  is  to  expand  upon  the  contents  of  these  sections, 
based  upon  the  e.xisting  hardware  which  differs  slightly  from  that  originally 
planned.  The  programmer  should  find  here  information  needed  to  write 
assembly  language  programs  for  the  purpose  of  communicating  between  the 
scan  converter  color  displays  and  the  analvzer  (Interdata  7/3J  minicomputer). 

The  first  section  consists  of  operating  instructions  for  the  LWCA 
control  panel  and  scan  conversion  processor,  while  subsequent  sections 
discuss  addressing  conventions  and  each  of  the  three  basic  types  of  data 
transfers:  Write  Display  M emory.  Read  Display  Memory,  and  Cursor  Data 
Entry.  Programming  examples  are  also  included. 

Operation  of  the  LWCA  Scan  Converter 

Scan  converter  operation  is  covered  in  Reference  2;  the  information 
presented  here  is  intended  to  serve  as  a supplement  and  co\ers  operation  with 
the  LWCA  Control  Panel  illustrated  in  Figure  1.  Except  for  the  ERASE 
DISPLAY  buttons,  all  of  the  switches  on  the  control  panel  also  ser\e  as 
indicators  controlled  by  their  state  and/or  the  DDI  (Display  Data  Interface) 
within  the  scan  converter.  An  exception  is  the  control  labeled  DAT.'X  SOURCE 
TAPE  which  functions  only  as  an  indicator  to  denote  the  fact  that  the  Precision 
Digital  Video  Integrator  has  been  set  to  accept  data  from  Mag  Tape  for  display 
on  the  scan  converter. 

The  l.WCA  TO  DISPL/XY  controls,  when  lit,  indicate  that  write 
and/or  read  data  transfers  are  enabled  in  the  hardware.  They  are  affected 
by  several  controls  on  the  Scan  Conversion  Processor  as  indicated  in  Table  1. 
The  state  of  these  controls  can  be  uniquely  determined  from  the  status  byte 
of  the  Display  Data  Port  which  has  been  assigned  device  number  X'8B'. 
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The  scan  converter  will  operate  normally  in  the  following  mode 
switch  positions:  PFl,  RHI,  CAPPI  and  B.  If  the  appropriate  I.WCA  TO 
DISPLAY  indicator  is  lit,  the  analyzer  can  read  or  write  into  the  display 
memories.  In  mode  switch  position  A,  the  necessary  conditions  for  I.WCA 
operation  are  set  up;  these  are; 

(1)  Scan  converter  in  RHI  mode, 

(<1)  LWCA  TO  DLSPLAY  RLAD  indicator  forced  ON 

(3)  LWCA  TO  DISPLAY  WRITK  indicator  forced  OK  if  the  converter 

memory  buss  is  available  (all  STORL  VIDEO  switches  Ol-'F). 
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The  scan  converter  ERASE  VIDEO  buttons  used  in  norn.al  operation 
do  not  erase  the  entire  screen;  the  contour  threshold  legend  area  is  left  un- 
changed. In  addition,  a mask  obscures  from  view  certain  areas  within  the 
ancillary  data  portion  of  the  screen.  These  areas  contain  coded  informtition 
available  to  the  analyzer  and  needed  by  the  contouring  hardware.  The  ER/\SE 
DISPLAY  buttons  on  the  l.WCA  CON  I'ROL  PANEL  not  only  erase  the  entire 
display,  but  also  iirhibit  the  mask  so  that  the  entire  screen  is  available  to 
display  information  from  the  processor.  The  mask  and  the  legend  are 
restored  when  the  operator  actuates  the  corresponding  STORE  TIIRESHOI.DS 
button  on  the  scan  conversion  processor. 

The  cursor  can  be  made  to  appear  in  any  display  by  depressing  the 
appropriate  CURSOR  ON/OFF  switch;  the  on  state  is  indicated  by  illumination 
of  tlie  switch.  The  cursor,  a blinking  single  point  on  the  dispLy,  can  be 
located  anywhere  on  the  screen  by  means  of  the  CURSOR  PO.SITION  trackball.^ 
The  cursor  changes  color  as  a function  of  its  surroundings  so  as  to  remain 
visible.  During  normal  scan  converter  operation,  the  mask  will  obscure  the 
cursor.  If  the  cursor  cannot  be  found,  the  following  property  may  be  useful: 
along  the  Top  and  Left  edges  of  the  display,  the  cursor  will  stop  even  if  the 
trackball  is  rotated  too  far.  At  the  bottom  edge,  the  cursor  disappears. 

When  moved  beyond  the  right  edge,  it  reappears  at  the  left  where  it  finally 
stops  about  an  inch  from  that  edge;  however,  if  the  SEND  DATA  button  were 
pressed  with  the  pursor  in  such  a position,  the  address  would  be  wrong. 

The  color / intensity  code  covered  by  the  cursor,  as  well  as  its 
coordinates,  can  be  entered  into  the  analyzer  by  pushing  the  appropriate 
SEND  DATA  button.  The  corresponding  cursor  must  be  switched-on  for 
this  action  to  be  recognized.  The  SEND  D.‘\ T A switch  will  light  when  de- 
pressed, if  the  DDI  control  logic  is  in  the  proper  state,  and  will  extinguish 
about  one-half  second  after  the  resulting  interrupt  has  been  serviced  by  the 
analyzer.^'  Pressing  the  INI  button  on  the  analyzer  console  should  always 
turn  off  anv  .SI-.N'D  DATA  indicators  which  are  lit  for  whatever  reason. 

Display  Conventions 

The  four  display  channels,  numbered  one  through  four,  contain 
independent  memories.  Each  memorv  is  organized  so  that  its  address 
corresponds  with  the  (X,  Y)  coordinates  within  a 248  by  320  point  matrix  as 
indicated  in  F'igure  2.  F.ach  point  can  take  on  one  of  sixteen  color/ intensity 
combinations  as  listed  in  that  figure.  (The  observed  colors  are  a function 
of  the  settings  of  an  array  of  switches  in  each  memory  interface  unit;  those 
colors  listed  correspond  to  the  settings  indicated  in  Figure  4-11  of  Ref.  2.  ) 
Note  that  color  15  has  a non-over  write  property;  once  this  code  occupies  a 
point,  the  only  way  the  color  code  at  that  point  can  be  changed  is  bv  erasure. 


^ If  any  .SEND  DAT/\  indicator  is  lit.  no  cursor  will  respond  to  the  trackball 
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The  ancillary  data  area  has  significance  only  in  normal  scan 
comerter  operation;  its  outline  is  indicated  in  Figure  Z,  while  the  details 
of  its  contents  appear  in  Figure  3.  Information  necessary  for  inter- 
pretation of  the  radar  video  data  portion  of  the  display  (scaling,  origin 
location,  time,  contour  thresholds,  and  antenna  angle)  is  obtainable  by 
reading  the  four -bit  codes  in  the  patches  indicated.  Fach  of  these  patches 
contains  the  same  four -bit  code  at  all  addresses  within  it.  Most  of  the 
code  patches  have  dimensions  of  5 .x  4 points  (the  same  as  the  color  patclies) 
e.xcept  for  the  origin  location  and  scaling  codes  which  are  only  5x1.  In 
either  case,  it  is  only  necessary  to  read  one  point  per  yjatch,  vinless  son’ie 
sort  of  error  correcting  scheme  is  implemented  to  make  use  of  the  redun- 
dancy. 


Points  written  as  color  15  by  a normally  operating  scan  converter 
(not  through  the  display  data  port)  within  the  ancillary  data  area  do  not  have 
the  non-o\erwrite  property.  Any  address  in  the  ancillary  data  area  which  is 
not  occupied  by  a 4 .x  5 patch  or  an  8 .x  5 character  can  be  used  for  storage 
of  a 4 -bit  word  (e.  g.  , to  "mark"  a stored  video  image)  except  for  the  8 .x  5 
area  under  each  color  patch.  Only  the  characters  and  color  patches  are 
displayed;  everything  else  in  the  ancillary  data  area  is  masked.  Again,  the 
entire  display  area  is  erased  (changed  to  color  zero)  and  the  mask  is 
inhibited  when  an  KRASF  DISPLAY  button  is  pushed,  the  entire  area  is  now 
a\ailable  to  accept  data  from  the  analyzer. 

General  Comments  on  the  Display  Data  Port 

The  hardware  which  comprises  the  display  data  port  controller 
consists  of  two  parts;  an  Interdata  Universal  Logic  Interface  (ULI)  and  a 
Raytheon -designed  Display  Data  Interface  (DDI).  The  ULI  responds  to 
device  address  X'8B'  and  contains  interrupt  and  byte/halfword  logic  controlled 
by  bits  0,  1 and  Z of  the  command  byte  (see  Figure  4,  note  Z).  Bit  Z should 
always  be  zero  since  the  display  data  port  operates  only  in  the  byte  mode. 

Bits  0 and  1 affect  interrupts  in  the  following  way;  0 1 -inter rupts  enabled; 

1 0 -i nter r upts  disabled  but  queued;  1 1 -inter rupts  disarmed  (neither  accepted 
nor  queued);  00 -previous  interrupt  state  unchanged.  The  ULI  does  not  affect 
any  bits  in  the  status  byte. 

The  DDI  contains  control  logic  which  is  described  by  the  state 
diagram  in  Figure  4.  Much  of  the  notation  here  will  not  be  of  concern  to  the 
programmer.  It  is  sufficient  to  note  that  state  transitions  are  typically 
caused  by  execution  of  the  7/3Z  I/C  instruction  listed  before  the  comment 
under  each  transition,  or  by  a hardware-generated  interrupt.  Operation  of 
the  DDI  control  logic  depends  on  the  state  of  bits  4,  5 and  6 of  the  command 
byte  as  tabulated  at  the  lower  right  of  Figure  4.  Also  located  there  is  a 
definition  of  the  status  byte,  of  which  bits  3 through  7 are  used. 
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Write  Display  Memory 

Three  distinct  types  of  write  operations  which  might  be  useful  in 
various  situations  are  supported  in  the  DDI  control  logic.  Controlled  by 
bits  4,  5 and  6 of  the  command  byte  (Figure  4),  they  include: 

(1)  0 0 0 - Write  single  point  or  multiple  points  the  same  color.  The 

first  write  instruction  transfers  and  color  code,  while  succeeding 

pairs  of  instructions  transfer  (X.\,  Y^).  The  notation  used  here  is  explained 
in  Figure  Z;  and  the  relationship  to  the  Interdaia  bit  nvunbers  can  be  determined 
from  Table  Z.“'  This  type  of  transfer  might  be  useful  where  many  points  of  the 
same  color  are  to  be  plotted  and  it  is  not  convenient  to  re-write  Sa.  XaM  and 
the  color  code  for  each  point.  An  example  is  listed  in  Table  Z.  After  the 
initial  write  instruction,  the  following  pairs  correspond  to  halfwords  so  that  a 
halfword  table  containing  (Xa,  Y a)  Nalues  could  be  easily  accessed  sequentially 
using  a write  block  instruction. 

(Z)  0 0 1 - Write  single  point  or  multiple  points  different  colors.  This 

sequence  operates  as  the  one  described  above,  except  that  after  the  Ya  transfer, 
the  next  instruction  transfers  another  nvmiber  for  Sa,  Xam  and  color. 

(3)  0 1 0 - Write  multiple  points,  fullword  boundaries.  This  sequence 

operates  as  the  one  described  above,  except  that  after  the  Ya  transfer,  the 
next  instruction  transfers  nothing  (sec  Figure  4,  state  W4),  while  the  one 
following  it  transfers  another  number  for  Sa,  XaMi  and  color.  This  type  of 
operation  is  intended  for  sequentially  writing  from  fullword  tables  where  each 
fullword  contains  Sa,  XajvT  COLOR,  Xa  and  Ya  for  one  point. 

The  Scan  Con\erter,  although  it  has  an  independent  memory  for  each 
display,  shares  a memory  address  buss  among  the  four  display  channels. 

When  one  or  more  STORE  VIDEO  switches  is  on,  this  buss  is  not  available  to 
the  display  data  port  and  the  write  display  memory  operation  is  disabled  in  the 
hardware.  It  is  also  disabled  for  certain  other  switch  settings  as  indicated  in 
Table  1.  Whenever  the  write  operation  is  disabled,  status  bit  5 is  zero.  Before 
a write  operation,  it  is  good  practice  to  check  status  to  determine  that  bit  5 is 
one,  although  nothing  will  happen  if  a write  is  attempted,  because  the  operation 
is  disabled  in  the  hardware.  Status  bit  4 should  be  checked  to  make  sure  it  is 
zero;  this  bit  indicates  that  a cursor  data  transfer  is  in  progress  and  that  the 
display  data  port  is  not  available. 


■ See  page  13  for  table. 
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Read  Display  Memory 

There  are  two  types  of  read  operations,  liependini’  on  whether  or 
not  tlae  scan  converter  memory  buss  is  a\aildble.  If  the  buss  is  available 
(all  STORK  VIDEO  switches  off;  status  bit  7-1),  then  a normal  read,  whicli 
operates  in  much  the  same  way  as  the  write  display  memory  transfer  des- 
cribed in  the  preceding  section,  can  be  executed.  Otlierwise,  the  process 
must  be  a slow  read,  which  involves  an  interrufit  serxice  routine.  Both 
types  of  read  operations  are  inhibited  if  the  read  indicator  is  not  lit  (status 
bit  6=0,  see  Table  1). 

Read  Display  Memory  --  Normal 


An  example  of  this  type  of  data  transfer  appears  in  Table  Z.  First, 
the  status  is  sensed  to  ensure  that  the  memory  buss  is  axailable,  the  read 
indicator  is  on,  and  that  no  cursor  data  entry  is  in  progress.  Next,  the 
proper  command  byte  is  output  to  device  X'8B'  and  Syy,  Xy\  and  Yyy 

are  transferred  just  as  for  the  write  operation.  At  this  point,  a delay  of  at 
least  si.xusec;  (for  e.xample,  four  BTCR  0,  0 (OZOO)  instructions)  must  be 
executed  so  that  the  hardware  is  sure  to  have  the  required  data  ready. 
Lesser  delays  might  work  but  have  not  been  tried.  Ne.xt,  a read  instruction 
is  e.xecuted;  the  4 -bit  color  code  appears  in  the  four  least  significant  bits  of 
the  second  operand.  Finally,  a command  byte  can  be  output  to  leave  the 
control  logic  in  state  I. 

Read  Display  Memory  --  Slow 

If,  in  the  preceding  section,  status  bit  7 had  been  found  to  be  zero, 
then  a slow-read  operation  must  be  used.  An  e.xample  is  found  in  Table  Z. 
Steps  0 through  6 are  the  same  as  for  a normal  read,  except  that  interrupts 
are  enabled.  The  control  logic,  after  step  6,  ends  up  in  state  SR4  (see 
F'igure  4)  where  it  waits  for  an  interrufit.  This  wait  could  last  as  long  as 
1 6 milliseconds  and  ends  when  the  DDI  has  obtained  data.  Other  processing 
can  be  executed  during  this  wait  interval.  When  the  interrupt  occurs,  a 
simple  interrupt  serxice  routine  consisting  of  steps  8 through  10  of  the 
example  in  Table  Z completes  the  operation. 


Details  on  interrupt  processing  can  be  found  in  Interdata  Docvurients; 

Model  7/3Z  Reference  Manual,  Pub.  No.  Z9  - 399R0Z,  Section  Z.  4 
3Z-Bit  Series  " " Pub.  No.  Z9  - 365R01,  Chapter  7. 
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As  does  the  slow-read  operation,  the  cursor  data  entry  makes 
use  of  an  interrupt  service  routine  and  a data  acquisition  method  which 
does  not  require  the  scan  con\erter  memory  buss.  There  is,  however, 
no  long  delay  because,  following  the  pressing  of  a SEND  DATA  button,  no 
interrupt  is  generated  until  after  all  required  data  has  been  obtained.  The 
cursor  data  entry  requires  that  the  DDI  control  logic  be  in  state  I and  that 
the  U LI  has  interrupts  enabled;  hence,  step  0 of  the  example  in  Table  Z. 

The  remainder  of  this  example  is  an  interrupt  service  routine  which  checks 
the  status  to  see  that  the  interrupt  was  caused  by  a cursor  data  entry, 
outputs  a command  byte  to  disarm  further  interrupts,  then  transfers  Sq, 
XcM>  -'*^C  second  operand  locations  of  the  three  read  in- 

structions. Finally,  the  control  logic  is  returned  to  state  I,  interrupts  are 
again  enabled,  and  the  original  program  status  word  is  restored. 

Programming  E.xamples 

The  SCPLF  subroutine  listed  in  Table  3 was  used  in  the  Liquid 
Water  Content  display  subroutine  to  take  care  of  getting  the  ninth  bit  of  X 
in  the  right  place  and  to  e.xecute  the  necessary  lO  instructions  for  writing 
one  point.  The  inputs  were  left  in  registers  and  the  subroutine  was  called 
using  BAL  F,  SCPLT.  Because  no  other  data  transfer  modes  were  being 
used  in  this  application,  the  command  byte  was  programmed  to  alwavs  leave 
the  ULIwith  interrupts  disarmed.  SCPLT  is  called  many  times  during  the 
main  program;  it  always  leaves  the  control  logic  in  state  I.  But  in  order  to 
ensure  that  the  very  first  point  is  plotted,  the  following  instructions  should 
be  executed  before  SCPLT  is  called  for  the  first  time: 

Ull  B,  X'8B' 

OC  B,  DDICMDZ 

Thus  forcing  the  control  logic  to  state  I. 

Another  way  to  structure  SCPLT  would  put  what  is  now  line  7Z 
(Table  3)  after  line  58.  thus  .SCPLT  would  not  leave  the  control  logic  in 
state  I,  but  would  force  it  there  first  each  time  it  is  called.  A third  method 
would  involve  forcing  the  control  logic  to  state  I only  once,  then  not  using 
any  OC  instructions  at  all  in  SCPLT  itself.  This  method  is  the  simplest 
and  fastest,  but  depends  on  nothing  disturbing  the  control  logic  between  calls 
of  SCPLT,  where  it  would  be  left  in  state  RW  1 (Figure  4). 
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Table  4 lists  a program  to  copy  one  display  to  another.  It  was 
written  directly  in  machine  language  as  a diagnostic  to  test  the  hardware, 
which  it  does  very  well  since  it  accesses  all  display  memory  locations  and 
exercises  the  read  circuitry  in  the  source  display  and  the  write  circuitry 
in  the  output  display.  A good  test  of  the  hardware  would  consist  of  the 
following : 

1)  Store  a test  pattern  or  radar  data  image  which  contains  all 
16  colors  in  display  1;  erase  displays  2,  3 and  4. 

2)  Put  cl  = 0 0 0 0 and  P = 0 0 2 0 into  the  program  and  run. 

Displays  1 and  2 should  now  be  identical. 

3)  Put  Q = 0 0 2 0 and  P = 0 0 4 0 into  the  program  and  run. 

Displays  1,  2 and  3 should  now  be  identical. 

4)  Put  0 = 0 0 4 0 and  P = 0 0 6 0 into  the  program  and  run. 

All  displays  should  be  identical. 

5)  Erase  display  1. 

6)  Put  Q = 0 0 6 0 and  P - 0 0 0 0 into  the  program  and  run.  All 
displays  should  again  bo  identical. 

Table  4 is  shown  set  up  for  a normal  read;  to  exercise  the  slow 
read,  follow  the  directions  at  the  end  of  the  table.  Execution  of  the  copy 
program  takes  about  three  seconds  in  the  normal  read  mode,  and  five 
seconds  in  the  slow  read  mode.  If  the  full  16  milliseconds  delay  were 
incurred  at  every  point  in  the  slow  read  mode,  the  program  would  require 
over  2 1 minutes  for  e.xecution  The  reason  it  only  takes  5 seconds  lies  in 
the  format  adopted  for  scanning  in  the  copy  program.  Examination  of 
Table  4 will  reveal  that  the  copy  process  is  basically  accomplished  by 
reading  one  point  from  the  source  display,  writing  that  data  into  the  same 
address  in  the  output  display,  incrementing  by  one  to  the  next  Y address,  then 
repeating.  When  Y reaches  248,  X is  incremented  by  one  and  Y goes  back  to 
zero.  The  fact  that  Y changes  more  rapdily  then  X is  the  key  to  the  reason 
for  the  unexpectedly  fast  performance  in  the  slow  read  mode.  The  average 
delay  is  only  about  67  microseconds  because  of  the  way  in  which  the  copy- 
progrcim  scan  interacts  with  the  display  raster-scan. 
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(X,  Y);  0 ^ X S 319;  0 « Y 

«247 

L 

(X,  8) ' 

Radar  Video 

Ancillary 

(254,  Y)  — 

Display  Area 

Data  Area 

(X,  179)-' 

(249,  Y)— 

(0,247)  (319.247) 


Color  Code  Relative  Video  Voltage 

(4-Bits)  Red  Green  Blue  Observed  Color 


Represented  by  8-bit  codes 
Display  coord.  - X^, 

Cursor  coord.  - 


Represented  by  1-bit  code: 
Display  coord.  - X^^ 

Cursor  coord.  - 

(The  most  significant 
bit  of  X) 

Represented  by  2-bit  codes: 

Display  select  - 

Cursor  select  - 

(Selects  one  of  the  four 
channels;  00=  display  If], 
01  = display  »2,  etc.  ) 
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Table  -1. 

P rug  null 

to  Cojiy  t)nc  Display  to  Another 

6000 

C8A0 

LHl 

A, ft 

source  display;  1254 

2 

ft  ; 0000  0020  0040  0060 

4 

C8B0 

LHI 

B,P 

output  display;  1234 

6 

P ; 0000  0020  0040  00 6o 

8 

C880 

LHI 

e,8B 

device  code  in  reg.  8. 

A 

00  8B 

C 

C890 

LHI 

9,C8 

cmd  byte  for  read  in  reg.  9* 

£ 

00C8 

6010 

C850 

LHI 

5,C2 

and  byte  for  write  in  reg.  5. 

2 

0002 

4 

24D0 

LIS 

D,0 

zero  reg.  D, 

6 

24E0 

LIS 

E,0 

M It  11  g 

8 

9E89 

OCri 

8,9 

output  cmd  byte  for  read. 

A 

9A8A 

8, A 

write  S ,X 

C 

9A8D 

\fOA 

6,D 

A 4UM  ^ ^ 

£ 

9A8E 

'WDR 

8,E 

" "^A 

6C20 

0200 

3TGR 

0,0 

delay  ^ "Vo  oe) 

2 

0200 

bTCK 

0,0 

M 

4 

0200 

ETCH 

0,0 

H 

0 

o;.oo 

ETCH 

0,0 

»( 

6 

9B6C 

HDH 

6,0 

read  data  into  reg.  C. 

A 

C4C0 

Ntil 

C,P 

mask  all  but  the  4 lab  of  reg.  C. 

C 

000  F 

£< 

9E85 

OCR 

8,5 

output  and  byte  for  write. 

b050 

OACE 

ALi 

C ,E 

get  the  output  display  code  in  reg.  C. 

2 

9A6C 

OR 

e,c 

write  Color 

4 

9A8D 

rfDii 

8,D 

II  ill  4 „ 

■^A 

6 

9A8E 

aER 

e,E 

" ’^A 

8 

2(E1 

AIS 

E,1 

increment  by  1. 

A 

C9E0 

CKI 

00 

U, 

com^-are  Ya  with  248 

C 

OOFS 

V 

4340 

BNP 

6018 

if  Ya  ^ 248,  go  to  £016. 

6040 

4000 

2 

bCl8 

» 

26D1 

AIS 

2,1 

increment  Xa  by  1 . 

6 

24E0 

LIS 

E,0 

zero  Ya. 

8 

C9D0 

CHI 

D,1C0 

compare  Xa  with  256. 

A 

0100 

c 

4210 

RI'. 

6018 

if  X^<  256,  go  to  6018. 

£ 

4000 

6050 

6016 

2 

4350 

ENE 

60 14 

if  Xa^256,  go  to  6Cb4. 

4 

4000 

6 

6064 

8 

C9D0 

CHI 

D,13F 

compare  Xa  with  519* 

A 

01 3F 

C 

4340 

BNF 

6018 

if  Xa^:  519,  go  to  6018. 

E 

4000 

6060 

6018 

f . 

mA 
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Table  4 

. (Continued) 

6062 

2200 

BFBS 

0,0 

branch 

unc.  to 

self. 

4 

6 

CAAO 

0010 

AHI 

A, 10 

add  1 6 

to  reg. 

A (make  ), 

8 

A 

CABO 

0010 

AHI 

B,10 

M It 

H II 

B " '* 

C 

E 

6070 

4500 

4000 

6018 

B 

ooie 

go  to 

601  8. 

To  do  the  same  task  iising  the  slow  read  mode,  change: 

600E  004A  cmd  hyte  for  slow  read,  interrupts  enabled. 

6020  2200  BFBS  0,0  branch  unc.  to  self. 

and  include  the  following  interrupt  service  routine: 

6100  2612  AIS  1 ,2  increment  reg.  1 , the  loc  part  of  the  PSW, 

by  one  halfword  to  bypass  the  2200  at  6020. 

6102  1 800  IKiWa  0,0  restore  the  PSW 


Hun  with  immediate  interrupts  enabled,  in  reg.  set  0.  (iSW=4000') 
In  the  interrupt  service  pointer  table,  at  Du+2x6B,  put  the 
starting  address  of  the  interrupt  service  routine: 

0E16  6100 


Note:  This  program  was  written  diectly  in  machine  langaitge;  it  wtis  never 
assembled  by  CAL.  The  assembler  notation  included  here  is  incorrect  for 
CAL  in  that  all  numbers  listed  are  in  hex.  In  CAL,  such  numbers  must  be 
represented  as  X'l.'IO.'N'  cr  Y ' lOaiXN'M.C ' , except  for  0-9* 
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The  LWCA  minicomputer  to  radar  interface  consists  of  two  software 
groups,  the  video  interface  and  the  ancillary  data  interface.  These  are 
operative  both  for  data  input  (recording)  and  output  (playback). 

The  interface  programs  have  been  implemented  as  general  purpose 
"I/O  drivers"  to  allow  other  LWC  system  modifications  without  any  changes 
to  these  modules.  There  are  4 entry  points  to  the  subroutines,  one  each 
for  video  input,  video  output,  ancillary  input  and  ancillary  output.  The 
calling  sequence  for  any  of  the  4 is  as  follows:- 

BAL  RX,  SUBRTNE  CXDTO  DESIRED  SUBROUTINE 

DC  Z(ERREXIT)  TRANSFER  ADDRESS  IF  ERROR 

OCCURED 

The  desired  entry  point  name  is  used  in  place  of  SUBRTNE,  and  the  desired 
return  address  register  is  used  for  RX.  If  no  error  occurs  in  the  execution 
of  the  subroutine,  control  is  returned  to  the  location  following  the  "DC" 
statement  w-hen  the  subroutine  operation  is  complete.  If  an  error  is  detected 
control  is  transferred  to  the  location  whose  label  is  substituted  for  ERREXIT, 
and  execution  continues  from  that  point.  When  an  error  occurs,  register 
10  contains  an  error  flag,  to  be  defined  below. 

The  individual  entry  points  are  set  up  as  follows: 

Video  Input 

Entry  name:  VIDINP 

Number  of  bytes  to  transfer:  in  location  RNUMBIN 

Location  for  start  of  transfer:  VBUF  + 32 

Note:  Actual  transfer  starts  at  VBUF  + 20,  to  allow  for  synchronization 

with  interface  hardware.  Actual  useful  daU  starts  at  VBUF  + 32  when  transfer 
is  complete. 

Error  List: 

0 No  e r ro  r 

1 No  selector  channel  start  - 6 second  timeout 

2 No  selector  channel  stop  - 2 millisecond  timout 
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Selector  channel  status  non-zero  when  transfer  complete  - 
probable  memory  error 

4 Full  transfer  not  completed  - abnormal  end 

5 Selector  channel  being  busy  to  start  - 6 second  timeout 

Subroutines  used: 

Transfers  control  to  VIDIST 

Sets  up  read  command  and  transfers  control  to  VIDI 
General  purpose  video  I/O  routine 
Universal  logic  interface  interrupt  routine  ("dump" 
pulse  interrupt) 

Selector  channel  end  of  operation  interrupt  routine 
(end  of  transfer  interrupt) 

Special  requirements:  The  actual  Input  or  Output  transfer  from  or  to  the 
radar  hardware  is  held  up  until  the  "dump"  pulse  is  received,  then  proceeds 
at  near  memory  speed.  ]f  no  dump  signal  is  received  for  approximately 
6 seconds  after  the  request  for  transfer,  an  error  exit  is  made. 

Video  Output 

Entry  name:  VIDOUT 

Number  of  bytes  to  transfer:  see  VIDINP 

Location  for  start  of  transfer:  see  \TDINP 
Error  List:  see  VIDINP 
Subroutines  used: 

VIDOUT:  Transfer  control  to  VIDOST 

VIDOST:  Sets  up  write  command  and  transfers  control  to  VID^ 

VIDl  to  ISELCH:  Same  as  VIDINP 

Special  requirementsr  same  as  VIDINP 

Ancillary  Data  Input 
Entry  name:  ANCINP 
Number  of  bytes  to  transfer:  8 
Location  of  start  of  transfer:  V^BUF 

Error  List: 

0 No  error 

1 No  data  available  interrupt  since  last  read/write 

2 Inputs  (or  outputs)  all  zero 


VIDINP 

VIDIST 

VIDl 

lULI 

ISELCH 
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ANCINP:  Transfer  control  to  ANCIST 

ANCIST:  General  Purpose  ancillary  data  I/O  setup  subroutine 

DELTAS:  Determine  direction  of  theta  and  phi  motion  (positive  or  negative) 

DELTAFIX;  Smooth  average  theta  and  phi  differences,  pass  to  pass 
R.tLI2I:  TTLI  interrupt  processor  for  reading  ancillary  data  via  the 

multiplexor  channel  in  response  to  the  dump  pulse. 

Special  requirements:  Exits  on  error  if  no  dump  pulse  is  received  between 

calls  (no  new  data  available)  or  if  the  data  is  all  zeros,  an  illegal  condition. 
The  theta  and  phi  directions  on  each  pass  are  differenced  from  the  previous 
pass  and  averaged  to  provide  an  indication  of  "positive  "or  "negative" 
position  change,  with  some  hysteresis  to  allow  for  indicator  wobble  . These 
flags  are  used  to  set  direction  bits  in  0 and  0 in  VBUF  to  guarantee  the 
proper  "painting"  of  display  fields  when  consecutive  data  samples  are  not 
contiguous  on  the  display. 

Data  is  not  actually  transferred  when  the  routine  is  called,  but  rather 
by  the  interrupt  routine  (lULlZIor  lULI  20)  when  the  dump  pulse  occurs. 

Ancillary  Data  Output 

Entry  name:  ANCOUT 

Number  of  bytes  to  transfer:  8 

Location  of  start  of  transfer:  VBUF 
Error  List:  Same  as  ANCINP 
Subroutines : 

ANCOUT:  Transfer  control  to  ANCOST 

ANCOST:  Initializes  theta  and  phi  direction  flags  and  transfers  control 

to  lable  ANC3  of  ANCIST 

ANCIST;  General  purpose  ancillary  data  I/O  setup  subroutines, 

DELTAS  and 

DELTAFIX:  see  ANCINP 

IULI20:  ULI  interrupt  processor  for  writing  ancillary  data  into  the 

multiplexor  channel  in  response  to  the  dump  pulse 
Special  Requirements;  Same  as  ANCOUT 

All  of  the  above  subroutines  are  described  in  detail  in  a general  purpose 
program  design  language  (PDL)  in  Appendix  H. 
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The  following  references  should  be  made  to  aid  in  a complete  under- 
standing of  data  formats,  equipment  operation,  etc. 

Raytheon  Memo  JHT:75:27  Revised.  Liquid  Water  Content  (LWC)  Constants, 

Variables  and  TSible ^contains  VBUF  format. 

Interdata  Manual  29-399  7/32  User's  Manual/ Conta ins  extended 
selector  channel  interface  description. 

Interdata  Manual  29-311  Universal  Logic  Interface  Instruction  Manual/ 

ULI  interface  description. 
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appendix  h 

LWC  Radar  Input/Output  Drivers 
Outlined  in  Program  Design  Language 


*VID15T 

VLDIST 


VIDl 


(VID9)' 


(VID8) 


(VID3) 


(VID2) 


(VID4) 

(VID7) 


VIDEO  INPUT  SUBROUTINE 
Save  registers  to  be  used  in  subroutine 
Get  read  Command 
Go  to  VIDl 

End  Subroutine  VIDIST 
Save  Command 

Save  Last  Byte  address  for  transfer 
Save  ULI  Interrupt  routine  address  into  trap 
Set  up  registeis  for  status  check 
Do  while  selector  channel  busy 

^6  second  timer  has  expired 

then  flag  selector  channel  hung  (5) 

Go  to  VIDERR 
Else  decrement  TIMBR 

Save  selector  channel  interrupt  routine  address  in  tra 
Output  a selector  channel  stop 
Flag  no  interrupts  recehred 
Enable  synch  interrupt 

Do  until  synch  interrupt  received 

6 8 econd  timer  has  expired 

then  flag  no  synch  interrupt  (1) 

Go  to  VIDERR 
Else  decrement  timer 
until  selector  channel  transfer  is  complete 

millisecond  timer  has  expired 

then  flag  no  selector  channel  stop  (2) 

Go  to  VIDERR 
Else  decrement  timer 

J^selector  channel  status  is  not  "successful  transfer" 
Then  flag  error  in  transfer  (3) 

Go  to  VIDERR 

^last  byte  address  not  desired  L.  B.  A. 

Then  flag  error  - Abnormal  end  (4) 


1. 


Labels  which  are  not  required  in  PDL  are  enclosed  in  parenthe 


VIDERR 


♦lULI 

lULI 


♦BELCH 

BELCH 

♦ANCIST 

ANCBT 

ANC3 

(ANC4) 

(ANC2) 

ANCERR 


Clear  Selector  Channel 

Restore  the  original  register  values 

If  transfer  had  no  errors 

Then  return  (to  call  point) 

Else  return  (to  error  exit  point) 

End  subroutine  VIDOST 

Universal  Logic  Interface  (ULI)  interrupt  Processor 

Turn  off  ULI  interrupt 
Start  Selector  channel  I/O 
Flag  synch  interrupt  received 
Return 

_End  subroutine  lULI 

Selector  Channel  Interrupt  Processor 

Save  selector  channel  status  as  done  flag  return 

End  subroutine  BELCH 

Ancillary  Data  Input 

Save  registers  to  be  used  in  this  subroutine 
Call  Deltas  (update  direction  flags  in  0 and  <2 ) 

Get  ULI  2 input  interrupt  processing  routine  address 

Save  in  ULI  2 interrupt  trap 
Turn  on  ULI  2 Interrupt 
If  first  pass  through  routine 

Then  Clear  ancillary  data  buffer  area 

Else  if  interrupt  has  not  been  received 
Then  flag  no  interrupt  received  (1) 

Go  to  ANCERR 
If  data  is  all  zero 

_^li£n_  flag  bad  data  (2) 

Clear  interrupt  received  flag 
Restore  Registers  to  origi  >1  values 
^ no  errors  occurred 

Then  return  (to  call  point) 

Else  return  (to  error  exit  point) 

Ejid  subroutine  ANCBT 
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♦ANCOST  Ancillary  Data  Output  Routine 
ANCOST  Save  registers  to  be  used  in  this  routine 
If  not  first  pass 

The_n  Call  DELTAS  (to  update  O and  0 Positions) 
If  Theta  direction  is  positive 

Then  flag  "positive"  in  ANCMSK+2 
Else  flag  "negative"  in  ANCMSK+2 
Get  address  of  output  ULI  interrupt  routine 
Go  to  ANC3  (in  routine  ANCIST) 

Else  set  Q and  0 directions  to  unknown 

Save  current  Q and  0 as  "last"  values 
Go  to  ANC6 

End  subroutine  ANCOST 

♦ Deltas  - Subroutine  to  find  the  new  Q and  0 positions 

DELTAS  Re  = current  O - last  0 (THETAL) 

Call  DELTAFIX  ( to  compute  average  difference) 
DTHEA  = Average  difference  in  O 
If  O direction  = negative 
Then  THEDIR  = - 1 
Else  THEDIR  = 1 
Re  = Current  0 - Last  0 (PHIL) 

_^a_l^  DELT AF IX  (to  compute  0 average) 

DPHIA  = Average  Difference  In  0 
If  0 direction  = negative 
Then  PHIDIR  = -1 
Else  PHUDIR  = 1 

THETAL  = Current  0;  PHIL  = Current  0 
Return 

End  subroutine  DELTAS 

♦ DELTAFIX  - Subroutine  to  average  differences 

DELTAFIX  RB  = Sign  flag  for  average  (O  or  0 ) 

If  sign  = minus 

Then  compliment  average 
RF  = average/4 
If  RF  Not  = 0 

Then  if  sign  was  negative 

Then  compliment  RF 


DFIX2 


RE  = RE  - RF  (subtract  average  from  current  value) 
Return 


*IULI2I 

1ULI2I 

lULl 

*rULl  20 
IULI20 


End  Subroutine  DELTt^FIX 

ULI  Input  interrupt  routine 

VBUF  = Data  from  Ancillary  Data  Port 
K_servo  angle  90°  (VBUF  + 6) 

Then  set  angle  = 0° 

ULINF  = '8000'  (flag  interrupt  occurred) 
Return 

_Ejif^  sub  routine  1UL12I 

ULI  2 output  interrupt  routine 
Set  (3  Direction  bit  from  ANCMSK 
Output  data  from  VBUF 
Clear  extra  Bit  in  (3 
Go  to  lULl 

End  Subroutine  IULI2.0 
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ATTN:  Thomas  Meagher 

ATTN:  Jerry  I . Harper 

ATTN:  Flank  fl.  Shelton 


ARC,  Incorporated 

ATTN:  John  C.  Adams 

AVCO  Research  ^ Systems  Group 


John  E.  Stevens,  JlOO 
S.  Skemp,  JtOO 
William  G.  Reinecke,  KlOO 
George  «»ebt'r,  .’<^30 


Battelle  Memorial  Institute 
ATTN:  Richard  Castle 

ATTN:  Me»*wvn  R.  Vanderl’nd 

ATTN:  E.  Unger 

ATTN:  A.  Ptcifer 


The  Boeing  Company 

ATTN:  Brian  lempriere 

ATTN:  Robert  Holmes 


LFE  Corporation 

Lnvironmental  Analysis  Lab.  Division 
ATTN:  Marxel  Nathans 

Lockheed  Missiles  s Space  Co.,  Inc. 


ATTN: 

Robert 

Au 

ATTN: 

Donald 

A.  Price 

ATTN; 

Charles  M.  Lee 

AT  TN : 

Arthur 

Coll  ins.  Dept 

ATTN: 

Gerald 

T.  Chrusciel 

Lockheed  Missiles  and  Space  fo. 


ATTN; 

T.  R. 

Fortune 

Marietta  Aerospace 

ATTN: 

,^i  1 1 iam  A.  Gray, 

ATTN: 

James 

M.  Potts, 

ATTN: 

Joyce 

LeGarc,  MR 

ATTN: 

Laird 

hinnaird 

Brown  Engineering  Company,  Inc. 
ATTN:  Ronald  Patrick 


ATTN:  Ken  Kreyenhagen 

CALSPAN  Corporation 

ATTN:  M.  S.  Holden 

ATTN:  Romeo  A.  Delibcris 

Effects  Technology,  Inc. 

ATTN:  Robert  Wenqler 


McDonnell  Douglas  Corporation 


ATTN; 

R.  J.  Reck 

ATTN: 

Ken  Kratch 

Inc. 

ATTN: 

H.  Hui-wiCJ 

ATTN. 

L.  Cohen 

Mi'teorology  Research,  Inc. 

ATTN:  Wil 1 iam  D.  Green 

National  Academy  ot  Sciences 

ATTN:  National  Materials  Advisory  Board  for 

IXinald  G.  Groves 

Pacif ic*Sierra  Research  Corp. 

ATTN:  Gary  Lang 


DiSt-3 


t 


DEPARTMtNMiF  lEFEMbE  CONTRACTORS  (Continued) 


Physicdl  Sciences,  Inc. 

ATTN;  M.  S.  Finson 


Physics  International  Company 

ATT.N;  Ooc.  Cun.  for  James  Shea 


Prototype  Development  Associates,  Inc. 
AlT*i!  J.  t.  Dunn 
Ar".:  L.  HuUdck 

ATTN:  R.  Hogan 


R 0 Associates 

AITN; 

Albert  L.  Latter 

ATTN; 

"aymond  F,  Ross 

ATTN: 

r.  A.  field 

ATTN; 

Cyrus  P.  Knowles 

ATTN; 

Harold  L.  Brodo 

;? 

The  Rand  Corporation 

! 

ATTN: 

R.  Robert  Rapp 

kay tneon  Company 

j 

ATTN; 

Libra ry 

r 

Raytheon  Company 

i 

ATTN; 

J.  H.  Turner 

f 

ATTN; 

•V.  C.  Anderson 

r 

3cy  ATTN: 

Anthony  J.  Jagodn 

OLPARTMI  NT 

Pr  ^ONTrA'  tors  iConiinued) 

Science  Appl icdtions,  Tnc. 

ATTN 

: John  earner 

Science  AppI icdtions,  !nc. 

ATTN 

Lyle  Dunbar 

ATTN 

: Car)  Swain 

Srience  Af plications.  Incorporated 

• .TN 

Aillidm  Layson 

2cy  ATTN 

Jonn  CocFayne 

Southern  Research  Institute 

ATTN 

C.  D.  Pears 

Stanford  Research  Institute 

ATTN 

Donald  Curran 

AITN 

George  R.  Abrahanson 

ATTN 

Phil  ip  J.  Dolan 

Systems,  Science  and  Softiyare,  Inc. 

ATTN 

G.  A,  Gurtman 

TRW  Systems 

Group 

ATTN 

] . E.  Aider,  Rl-lOOii 

ATT.N 

Thomas  6.  williams 

ATTN 

w.  w.  Wood 

ATTT, 

R.  K.  Plebuch.  Rl-3076 

ATTN 

Peter  Brandt,  El-2006 

2cy  ATTN 

D.  H.  Baer,  R1-2I36 

TRW  Systems 

broup 

ATTN; 

William  Polich 

ATTN: 

Earl  w.  Allen,  520/141 

ATTN: 

L.  Berger 

ATTN: 

E.  V.  Wong,  527/712 

ATTN: 

V.  Blankenship 

